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Abstract 

Background: Genome-wide significant associations of schizophrenia with eight SNPs in the CNNM2, MIR137, 
PCGEM1, TRIM26, CSMD1 , MMP16, NT5C2 and CCDC68 genes have been identified in a recent mega-analysis of 
genome-wide association studies. To date, the role of these SNPs on gray matter (GM) volumes remains unclear. 

Methods: After performing quality control for minor-allele frequency > 5% using a JPT HapMap sample and our 
sample, a genotyping call rate > 95% and Hardy-Weinberg equilibrium testing (p > 0.01), five of eight SNPs were 
eligible for analysis. We used a comprehensive voxel-based morphometry (VBM) technique to investigate the effects 
of these five SNPs on GM volumes between major-allele homozygotes and minor-allele carriers in Japanese patients 
with schizophrenia (n= 173) and healthy subjects (n = 449). 

Results: The rs7914558 risk variant at CNNM2 was associated with voxel-based GM volumes in the bilateral inferior 
frontal gyri (right T=4.96, p = 0.0088, left T=4.66, p = 0.031). These peak voxels, which were affected by the variant, 
existed in the orbital region of the inferior frontal gyri. Individuals with the risk G/G genotype of rs7914558 had 
smaller GM volumes in the bilateral inferior frontal gyri than carriers of the non-risk A-allele. Although several effects 
of the genotype and the genotype-diagnosis interaction of other SNPs on GM volumes were observed in the 
exploratory VBM analyses, these effects did not remain after the FI/l/ZF-correction for multiple tests (p > 0.05). 

Conclusions: Our findings suggest that the genetic variant in the CNNM2 gene could be implicated in the 
pathogenesis of schizophrenia through the GM volumetric vulnerability of the orbital regions in the inferior frontal gyri. 

Keywords: Schizophrenia, Genome-wide association study, Voxel-based morphometry, Cyclin M2 {CNNM2), 
Inferior frontal gyrus 



Background 

Schizophrenia is a common and complex psychiatric dis- 
order with a lifetime risk of approximately 1%. This 
disorder has a strong genetic component; indeed, the esti- 
mated heritability is 81% [1]. Multiple genetic variants that 
have a small effect have been implicated in the pathogen- 
esis of schizophrenia [2]. A genome-wide association study 
(GWAS) of single-nucleotide polymorphisms (SNPs) that 
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accesses tens of thousands of DNA samples from patients 
and controls can be a powerful tool for identifying common 
risk factors for complex diseases, such as schizophrenia. To 
date, GWASs on schizophrenia have identified several 
genome-wide significant associated variants located in the 
zinc finger protein 804A (ZNF804A), neurogranin (NRGN), 
transcription factor 4 (TCF4) genes and a major histocom- 
patibility complex (MHC) region [3,4]. Subsequently, the 
influences of these SNPs in the genes on brain function and 
structure have been reported [5]. We have found that the 
genome-wide supported variant of the NRGN gene is asso- 
ciated with the brain morphology of the anterior cingulate 
cortex in patients with schizophrenia [6]. 
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Recently, a combining analysis of a mega-analysis of 
GWAS data from 17 separate studies (9,394 cases and 
12,462 controls) and the replication data (8,442 cases and 
21,397 controls) of European ancestry have found genome- 
wide significant associations of schizophrenia with eight 
SNPs [rs7914558 (cyclin M2; CNNM2), rsl625579 (micro- 
RNA 137, MIR137), rsl7662626 (PCGEM1, prostate-specific 
transcript, PCGEM1), rs2021722 (tripartite motif containing 
26, TRIM26), rsl0503253 (CUB and Sushi multiple do- 
mains 1, CSMD1), rs7004633 (matrix metallopeptidase 
16, MMP16), rslll91580 (5 '-nucleotidase, cytosolic II, 
NT5C2) and rsl2966547 (coiled-coil domain containing 
68, CCDC68)] from the five new (lp2L3, 2q32.3, 8p23.2, 
8q21.3 and 10q24.32-q24.33) and two previously reported 
(6p21.32-p22.1 and 18q2L2) loci [7]. In studies exploring 
brain activation during a sentence completion task, indi- 
viduals without the risk allele of the MIR13 7 genotype had 
significantly greater activation in the posterior right medial 
frontal gyrus than at-risk individuals [8]. To our know- 
ledge, however, no study has investigated the effects of 
these SNPs on gray matter (GM) volumes. 

Many attempts have been made to minimize the clinical 
and genetic heterogeneity in studies of schizophrenia. One 
strategy for gene discovery uses neurobiological quantita- 
tive traits (QT) as intermediate phenotypes that reflect the 
underlying genetic vulnerability better than diagnostic 
categorization, such as schizophrenia [9,10]. This strategy 
has the potential to reduce clinical and genetic heterogen- 
eity [11]. Structural GM volumes indicate substantial 
heritability rates ranging from moderate (40-70%) to high 
(70-95%) in the frontal and temporal brain regions 
[12,13]. Mata analyses of brain morphological studies in 
individuals with first-episode schizophrenia and neuro- 
leptic naive schizophrenia as well as chronic patients 
with schizophrenia revealed reduction of GM volume 
in frontal, striato-limbic and temporal regions were 
present in the early stage of schizophrenia and were un- 
related to the effects of neuroleptic treatment, chron- 
icity and duration of illness [14-16]. Some studies have 
shown that abnormalities in GM volumes are inter- 
mediate phenotypes that bridge the gap between the 
genotype and diagnostic categorization [11,17]. Charac- 
terizing the functional effects of novel and poorly 
understood genetic variants on the intermediate phe- 
notypes provides important insights into the neural 
mechanisms by which the variants increase the risk 
for schizophrenia [18]. Our research group has had a 
long-standing interest in the effects of genetic vari- 
ants (i.e., COMT, DISCI, PACAP, BDNF, APOE, AKT1 
and NRGN) on brain structure in psychiatric disor- 
ders [6,19-24]. In this study, we examined the impact 
of the genome-wide supported variants on the GM 
volumes of patients with schizophrenia and healthy sub- 
jects. Using a comprehensive voxel-based morphometry 



(VBM) technique, we tested the hypothesis that these risk 
variants would be associated with GM volumes. 

Material and methods 

Subjects 

VBM analyses were conducted on 173 patients with 
schizophrenia (59.0% males, 102 males and 71 females, 
mean age ± SD, range 36.0 ± 12.3 years) and 449 healthy 
subjects (47.7% males, 214 males and 235 females, mean 
age ± SD, range 35.4 ± 12.8 years). All subjects were bio- 
logically unrelated within the second-degree of relationship 
and were of Japanese descent [25,26]. The subjects were 
excluded if they had neurological or medical conditions 
that could potentially affect the central nervous system, 
such as atypical headache, head trauma with loss of con- 
sciousness, chronic lung disease, kidney disease, chronic 
hepatic disease, thyroid disease, active cancer, cerebrovas- 
cular disease, epilepsy, seizures, substance-related disorders 
or mental retardation. Patients were recruited from the 
Osaka University Hospital. Each patient with schizophrenia 
had been diagnosed by at least two trained psychiatrists ac- 
cording to the criteria from the Diagnostic and Statistical 
Manual of Mental Disorders, Fourth Edition (DSM-IV) 
based on the Structured Clinical Interview for DSM-IV 
(SCID). Controls were recruited through local advertise- 
ments at Osaka University. Psychiatrically, medically and 
neurologically, the healthy subjects were evaluated using 
the non-patient version of the SCID to exclude individuals 
who had current or past contact with psychiatric services 
or who had received psychiatric medication. Current symp- 
toms of schizophrenia were evaluated using the positive 
and negative syndrome scale (PANSS) [27]. Mean age and 
handedness did not differ significantly between the study 
group and the controls (/?>0.38), while the female ratio, 
years of education, estimated premorbid intelligence quo- 
tient (IQ) and total gray matter volumes were significantly 
lower in the patients with schizophrenia compared to the 
controls (p < 0.016) (Additional file 1: Table SI). All partici- 
pants provided written informed consent after the study 
procedures had been fully explained. This study was per- 
formed in accordance with the World Medical Associa- 
tions Declaration of Helsinki and was approved by the 
Research Ethical Committee of Osaka University. 

SNP selection and SNP genotyping 

We selected eight SNPs, rs7914558, rsl625579, rsl7662626, 
rs2021722, rsl0503253, rs7004633, rsl 1191580 and 
rsl2966547, from a previous mega analysis of GWASs 
[7]. Of these eight SNPs, rsl7662626 in the PCGEM1 gene 
(2q32.3) was not polymorphic in the samples obtained 
from the HapMap Japanese in Tokyo (JPT) project. Be- 
cause the other seven SNPs were common genetic variants 
in the HapMap JPT samples (minor allele frequency > 5%), 
we focused on these SNPs. Venous blood was collected 
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from the subjects, and genomic DNA was extracted from 
whole blood according to standard procedures. These 
SNPs were genotyped using the TaqMan 5'-exonucle- 
ase allelic discrimination assay (Assay ID: rs7914558: 

C_31978821_10, rsl625579: C 8946584_20, rs2021722: 

C_11690541_10, rsl0503253: C 1503810_20, rs7004633: 

C_29048976_10, rslll91580: C_31656012_10 and 

rsl2966547: C 152930_10, Applied Biosystems, Foster 

City, CA, USA), as previously described [19,22]. Detailed in- 
formation on the PCR conditions is available upon request. 
The genotyping call rates were 98.2% (rs7914558), 97.9% 
(rsl625579), 83.6% (rs2021722), 97.9% (rsl0503253), 97.9% 
(rs7004633), 99.8% (rslll91580) and 97.6% (rsl2966547). 
No deviation from the Hardy- Weinberg equilibrium (HWE) 
in the examined SNPs was detected in the patients or in con- 
trols (p > 0.01), with the exception of rs2021722. A significant 
deviation from HWE in the rs2021722 was found in both 
the patients (p = 1.02 x 10~ 17 ) and controls (p = 2.05 x 10~ 36 ) 
with a relative excess of CC homozygotes, TT homozygotes 
and undetermined subjects. According to the dbSNP data- 
base (National Center for Biotechnology Information), the 
SNP is shown as a tri-allelic variant with T/C/A. A number 
of genome-wide significant variants within MHC (6p21.32- 
p22.1), including rs2021722, have been identified [7]. How- 
ever, the MHC region has been excluded from further ana- 
lysis. Analyzing the region is difficult because of its high 
linkage disequilibrium (LD) and ethnic heterogeneity. Minor 
allele frequencies of rsl625579 were under 5% in our patients 
(3.2%) and controls (2.5%). Therefore, in this study, we ex- 
cluded these SNPs rs2021722 and rsl625579 from the VBM 
analyses. Genotype and allele distributions for each SNP in- 
cluded in the VBM analyses between the patients with 
schizophrenia and the controls are shown in Table 1. All risk 
alleles were defined based on the previous GWAS [7]: 
rs7914558 (major G-allele), rsl0503253 (minor A-allele), 
rs7004633 (minor G-allele), rslll91580 (major T-allele) and 
rsl2966547 (minor G-allele). To increase the statistical 
power and decrease type I errors, homozygotes and the het- 
erozygotes for the minor allele groups were combined and 
treated as minor-allele carriers. In this study, we contrasted 



GM volumes between minor allele carriers and major allele 
homozygotes. 

Magnetic resonance imaging procedure 

All magnetic resonance imaging (MRI) studies were 
performed on a 1.5 T GE Signa EXCITE system. A three- 
dimensional volumetric acquisition of a Tl -weighted gradi- 
ent echo sequence produced a gapless series of 124 sagittal 
sections using a spoiled gradient recalled acquisition in the 
steady state (SPGR) sequence (TE/TR, 4.2/12.6 ms; flip 
angle, 15°; acquisition matrix, 256 x 256; 1NEX, FOV, 24 x 
24 cm; slice thickness, 1.4 mm). Subjects with MRI abnor- 
malities, such as infarcts, hemorrhages or brain tumors, 
were screened out prior to including this study as part of 
routine clinical diagnosis and treatment. Therefore, there 
were no gross abnormalities in any of the subjects. Each 
image was visually examined to eliminate images with mo- 
tion or metal artifacts, and the anterior commissure- 
posterior commissure line was adjusted. The MRI images 
were processed using the VBM8 toolbox in Statistical Para- 
metric Mapping 8 (SPM8; http://www.fil.ion.ucl.ac.ul</spm/ 
software/spm8/) running on MATLAB R2013a (Math- 
Works, Natick, MA, USA) according to the VBM8-Toolbox 
Manual (http://dbm.neuro.uni-jena.de/vbm8/VBM8-Manual. 
pdf). The Tl images were normalized and segmented into 
GM, white matter (WM) and cerebrospinal fluid (CSF) using 
the VBM8 toolbox with defaults for the extended options. 
The modulated non-linear only (i.e., with no affine com- 
ponent) option was selected to create volumetric GM par- 
titions. Finally, the images were smoothed with an 8 -mm 
full-width, half-maximum isotropic Gaussian kernel. 

Statistical analyses 

In genetic association analysis, we performed power calcula- 
tions using the Power Calculator for Two-Stage Association 
Studies (http://www.sph.umich.edu/csg/abecasis/CaTS/) [28]. 
The power estimate was based on an allele frequency of 0.53 
at rs7914558, a prevalence of 0.01, an alpha level of 0.05, and 
assuming varying degrees of odds ratio using a multiplicative 
model. In brain morphological analyses, we performed power 



Table 1 Genotype and allele distributions for each SNP between the patients with schizophrenia and healthy subjects 









Risk 




Genotype frequencies 




Risk allele 


Allelic 


OR 








allele 


+/+ 


+/- 


-/- 


+/+ 


+/- 


-/- 


frequencies 


p value 


(95% CI) 


SNP IDs 


Gene 


Chr 


+/- 


SCZ (n = 


173) 


CON (n = 


449) 


SCZ 


CON 


1X2) 




rsl 0503253 


CSMD1 


8p23.2 


A/C 


0.06 


0.46 


0.49 


0.09 


0.46 


0.45 


0.29 


0.32 


0.22 (1.5) 


1.19 (0.90-1.56) 


rs7004633 


MMP16 


8q21.3 


G/A 


0.03 


0.37 


0.60 


0.06 


0.36 


0.57 


0.21 


0.25 


0.24 (1.4) 


1.20 (0.89-1.61) 


rs7914558 


CNNM2 


10q24.32 


G/A 


0.23 


0.53 


0.24 


0.29 


0.49 


0.22 


0.49 


0.53 


0.23 (1.4) 


1.16 (0.91-1.49) 


rsl 1 1 91 580 


NT5C2 


10q24.33 


T/C 


0.49 


0.42 


0.09 


0.50 


0.42 


0.07 


0.70 


0.72 


0.55 (0.4) 


1.09 (0.83-1.43) 


rs 12966547 


CCDC68 


18q21 


G/A 


0.15 


0.44 


0.41 


0.17 


0.44 


0.39 


0.37 


0.39 


0.54 (0.4) 


1.09 (0.84-1.41) 



Abbreviations: Chr Chromosome; SCZ patients with schizophrenia; CON healthy controls; +, risk allele; -, non-risk allele; OR odds ratio. 
For alleles, the first allele is the risk allele. All risk alleles are represented based on the previous GWAS [7]. 
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calculations using the G*Power Version 3.1.5 [29]. The 
power estimate was based on an alpha level of 0.05, a power 
of 0.80 and assuming varying degrees of effect size using 
t tests. Standardized effect sizes were calculated using 
Cohens d method (http:www.uccs.edu/faculty/lbecker). 

Statistical analyses of the demographic variables were 
performed using PASW Statistics 18.0 software (SPSS 
Japan Inc., Tokyo, Japan). Based on the assumption that 
most of demographic variables, such as age and education 
years, were not fitted to a normality distribution with the 
Kolmogorov-Smirnov test (p < 0.05), differences in clinical 
characteristics between patients and controls or between 
genotypes were analyzed using the non-parametric Mann- 
Whitney U-test for continuous variables, such as age and 
years of education, and tests for categorical variables, 
such as gender and handedness, as shown in Additional 
file 1: Table SI -S3. The presence of HWE was examined 
using the test for goodness-of-fit via SNPAlyze V5.1.1 
Pro software (DYNACOM, Yokohama, Japan). The allelic 
distributions of each SNP between patients and controls 
were analyzed using ^ tests with the SNPAlyze software. 
The significance levels for HWE and other statistical tests 
were set at two-tailed ^-values, j?<0.01 and ^><0.05, 
respectively. 

We performed a comprehensive exploratory whole brain 
search using the SPM8 statistical tools to examine the ef- 
fects of the diagnosis, the genotype and their interaction of 
each SNP on GM volume in patients with schizophrenia 
and healthy subjects. As two-way ANOVA can simultan- 
eously investigate these effects only in one model, these ef- 
fects on GM volume were statistically assessed using a full 
factorial model for a 2 x 2 ANOVA with diagnosis (cases 
and controls) and genotype status (major-allele homozy- 
gotes and minor-allele carriers) as independent variables in 
SPM8. Gray matter volumes are correlated with age, and 
gender and years of education differed significantly be- 
tween the patient and the control groups (Additional file 1: 
Table SI). Therefore, age, gender and years of education 
were included as covariates of no interest into the analyses 
to control for confounding variables. We contrasted GM 
volumes between the diagnostic groups (smaller volume 
region in patients with schizophrenia compared with 
healthy subjects), the genotype groups (smaller or larger 
volume region in minor-allele carriers relative to major- 
allele homozygotes) or the genotype-diagnosis inter- 
action. Non-sphericity estimation was used. We applied 
a voxel-level height threshold of p < 0.001 (uncorrected for 
multiple comparisons) and clusters of more than 100 con- 
tiguous voxels were considered for the exploratory VBM 
analyses. And then we applied family-wise error (FWE) cor- 
rection for multiple testing to avoid type I errors at the 
whole brain level. Eventually, the significance level was set 
at p<0.05 (FWE- corrected ). To obtain a cluster as large as 
possible, we extracted relative GM volumes from nominal 



clusters in bilateral inferiror frontal gyri at the lenient uncor- 
rected threshold of p < 0.001 and cluster sizes > 100. The 
extraction of these relative GM volumes were performed after 
including confounding factors such as age, sex and education 
years and modulated by total brain volumes in the VBM ana- 
lyses. Anatomic localization was performed according to both 
the MNI coordinates and Talairach coordinates, which were 
obtained from M. Bretts transformations (http://imaging. 
mrc-cbu.cam.ac.uk/imaging/MniTalairach) and presented as 
Talairach coordinates. 

Results 

Our study size of 173 cases and 449 controls had suffi- 
cient power (>80%) to detect a genetic effect at odds 
ratio of 1.43 or greater for rs79 14558 when the allele 
frequency was 0.53. Unfortunately, in our sample sizes, 
there was no allelic association with schizophrenia for 
any of the five SNPs [rs7914558 (CNNM2), rsl0503253 
(CSMD1), rs7004633 (MMP16), rslll91580 (NTSC2) 
and rsl2966547 (CCDC68)] (p > 0.22, Table 1). 

We investigated the effects of diagnosis (cases and con- 
trols), genotype (major-allele homozygotes and minor- 
allele carriers) and their interaction of five SNPs on GM 
volumes in a comprehensive exploratory VBM analysis. 
The effects of diagnosis between patients with schizophre- 
nia and healthy subjects were found in all analyses of the 
present study (p < 0.05). Patients with schizophrenia 
showed smaller GM volumes compared with healthy sub- 
jects primarily in the frontal and temporal lobes, including 
the bilateral inferior frontal gyri, which was consistent with 
previous studies [14,30]. We found significant effects for 
the risk-allele homozygotes of rs7914558 on decreased 
GM volume in the bilateral inferior frontal gyri (right, T - 
4.96, p = 0.0088; left, 7=4.66, = 0.031), as shown in 
Table 2 and the regions based on the hot color map in 
Figure 1. To compare the effects of the genotype in both 
the patients with schizophrenia and healthy subjects, we 
extracted the means and SD for relative GM volumes from 
nominal clusters in bilateral inferior frontal gyri and 
the extracted GM volumes were shown in Figure 2 and 
Additional file 1: Table S2. The risk G -allele homozygotes 
of the CNNM2 polymorphism had smaller GM volumes in 
the bilateral inferior frontal gyri compared to the non-risk 
allele carriers. The inferior frontal gyrus can be subdivided 
into three macroanatomical structures: the orbital (Brod- 
mann area: BA47), opercular (BA44) and triangular (BA45) 
parts. In the present study, the peak GM regions affected 
by the CNNM2 genotype existed in the orbital parts of the 
bilateral inferior frontal gyri. When the two genotype 
groups (G-allele homozygotes and A-allele carriers) were 
divided into three genotype groups (individuals with GIG 
genotype, G/A genotype and A/A genotype) and we per- 
formed an additional VBM analysis using a multiple regres- 
sion model, the number of risk G-allele was significantly 
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Table 2 Effects of the CNNM2 genotype and genotype-diagnosis interaction on GM volumes 



Brain regions 



R/L 



BA 



CS 



p values (peak) 
T FWE 



Talairach coordinates 



Non-risk minor allele carrier > Risk 


major allele homozygote 














Inferior frontal Gyrus 


R 


11/47 


1306 


4.96 


0.0088 


22 


31 


-20 


Inferior frontal Gyrus 


L 


47 


437 


4.66 


0.031 


-22 


18 


-22 


Middle frontal Gyrus 


L 


11 


667 


4.33 


0.11 


-24 


38 


-18 


Posterior cingulate 


L 


29 


248 


3.73 


0.61 


-7 


-48 


11 


Non-risk minor allele carrier < Risk 


major allele homozygote 














no suprathreshold clusters 


















Genotype-diagnosis interaction 


















Superior temporal Gyrus 


L 


22 


598 


4.23 


0.16 


-46 


-20 


1 



Abbreviations: R right; L left; BA Brodmann area; CS Cluster size; FWE family-wise error. 

All regions shown have nominal association at a voxel-level height threshold of uncorrected P < 0-001 and a minimum clusters extent of 100 voxels. Significant 
results {FWE- correctedl p < 0.05) are shown in bold face. 



related to smaller GM volumes of the right inferior frontal 
gyrus (T=4.67, p = 0.029) but not the left inferior frontal 
gyrus (T= 3.64, p = 0.70) in total subjects. 

As shown in Additional file 1: Table SI, patients with 
schizophrenia participated in this study were chronic and 
the symptoms were moderately stable. There was no dif- 
ference in demographic information, such as duration of 
illness or PANSS scores, between risk G -allele homozy- 
gotes and non-risk A-allele carriers of rs79 14558 (p > 0.20, 
Additional file 1: Table S3). However, to find whether there 
was potential clinical impact on our outcomes, we add- 
itionally investigated the genotype effect on GM volume 



including chlorpromazine equivalent of total antipsychotics 
(mg/day), duration of illness or PANSS scores as covariant 
for VBM analyses only in patients. The genotype effect on 
these regions did not change even after controlling for 
these factors, suggesting that there was no potential clinical 
impact on our outcomes. 

In the exploratory VBM analyses, we also found nominal 
effects of the risk-allele homozygotes of rs7914558 on de- 
creased GM volume in the left middle frontal gyrus and 
the left posterior cingulate and a nominal genotype- 
diagnosis interaction with the GM volume in the left su- 
perior temporal gyrus (uncorrected P< 0.001, Table 2). 




WW WW 




Figure 1 Effects of the rs7914558 polymorphism at the CNNM2 gene on GM volumes. There were effects of the risk-allele carriers of 
rs7914558 at CNNM2 on decreased GM regions (red areas shown on the hot color map). There was no effect of the genotype on increased GM 
regions (blue area shown on the winter color map). Each color map shows the t values corresponding to the color in the figure. 
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R 



(B) Left inferior frontal gyrus 



R 



S 0/70 



51 



£ 0.62 



O 0.58 



2 0.54 



■ G/G 

□ A carrier 



I 



Schizophrenia Control 



(C) Right inferior frontal gyrus 

0.64 



N 



O 0.52 



□ G/G 

□ A carrier 



Schizophrenia Control 



Figure 2 The impacts of the CNNM2 genotype on GM volume of the bilateral inferior frontal gyri. (A) Anatomical localizations are 
displayed on coronal, sagittal, and axial sections of a normal MRI spatially normalized into the Montreal Neurological Institute template (p <0.05). 
The most significant cluster of the genotype effect was in the right inferior frontal gyrus. The region is shown as a cross-hairline. The color bars 
show the t values corresponding to the color in the figure. (B, C) Each column shows the relative GM volumes extracted from a nominal cluster 
at uncorrected P < 0.00] and cluster size > 100 in the right inferior frontal gyrus (peak Talairach coordinates; 22, 31, -20) (B) and in the left inferior 
frontal gyrus (peak Talairach coordinates; -22, 18, -22) (C). Error bars represent the standard error. 



Additionally, we found several marginal effects of genotypes 
and genotype-diagnosis interactions of other SNPs on GM 
volumes in the exploratory analyses (uncorrected P< 0.001, 
Additional file 1: Table S4 and Figure S1-S4). However, 
these effects of genotypes and genotype-diagnosis interac- 
tions on GM regions did not survive after the FWE-correc- 
tion for multiple tests at the whole brain level (p > 0.05). 

Discussion 

To date, it remains unclear whether the genome-wide 
significant risk variants for schizophrenia in a mega- 
analysis of GWASs influenced GM volumes. This study 
is the first to identify the GM morphology associated 
with genome-wide risk variants using a comprehensive 
VBM technique. Of the five genetic variants investigated 
in this study, we found influences of the CNNM2 geno- 
type on the bilateral inferior frontal gyri at the whole 
brain level. GM volumes in the bilateral inferior frontal 
gyri, particularly the orbital region, in the risk G-allele 
homozygotes of CNNM2 polymorphism were smaller 
than those observed in the non-risk allele carriers. 

In genetic association analysis, our sample size had suffi- 
cient power (>80%) to detect a genetic effect with odds ra- 
tio of 1.43 or greater for rs79 14558. However, previous 
large GWAS has reported the genetic effect with low odds 
ratio of 1.10 for the SNP [7]. To detect such a small gen- 
etic effect, our sample size had insufficient power (12%), 



and a large sample size of at least 3400 patients and 3400 
controls is needed. On the other hand, in brain morpho- 
logical analyses, our sample size had sufficient power 
(>80%) to detect a genotype effect on GM volumes at 
medium effect size (Cohen s d) of 0.25 or greater. The ob- 
served effect sizes on the inferior frontal gyri were 
medium to large (0.29-1.00). These findings suggest that 
the genome-wide supported variant of schizophrenia had 
larger effect on GM volumes of the inferior frontal gyri 
than diagnostic status, and support that GM volumes ab- 
normalities were prominent intermediate phenotypes 
bridging the gap between a susceptibility genetic variant 
and diagnostic categorization. 

Rs7914558 is located in intronl of the CNNM2 gene 
(also known as ACDP2) on chromosome 10q24.32. The 
CNNM2 gene spans 160.3 kb of genomic DNA and con- 
tains eight exons. This gene belongs to a member of the 
ancient conserved domain-containing protein family be- 
cause the protein shares a domain conserved in a large 
number of species ranging from bacteria to human [31]. 
Members of this protein family contain a sequence motif 
that is present in the cyclin box, a cyclic nucleotide- 
monophosphate (cNMP) -binding domain. The CNNM2 
gene has a ubiquitous expression pattern in humans [31]. 
In particular, the level of CNNM2 expression in the 
brain is moderate to high (http://www.ebi.ac.uk/gxa/ 
experiment/E-MTAB-37/ENSG00000148842). 
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However, whether the expression level of this gene in the 
brains of patients with schizophrenia is lower or higher 
than that in healthy subjects is unknown. The encoded 
protein CNNM2 plays an important role in magnesium 
homeostasis by modulating Mg2+ concentration. The 
CNNM2 mRNA is upregulated when there is a defi- 
ciency of magnesium in the brain [32]. CNNM2 medi- 
ated Mg2 + -sensitive Na + currents were blocked by 
increased extracellular Mg2+ concentrations [33]. We 
assessed the effect of the rs79 14558 genotype on 
CNNM2 expression using bioinformatics data (http:// 
www.sanger.ac.uk/resources/software/genevar/ [34] ) to 
examine whether the rs7914558 genotype might be an 
expression quantitative trait loci (eQTL). In silico ana- 
lysis showed that the CNNM2 gene expression of the 
high-risk G genotype of rs79 14558 was significantly 
lower than that of the non-risk genotype in the com- 
bined lymphoblast-derived HapMap CEU and YRI sam- 
ples (r= -.23, t = -2.59, p = 0.011). The low expression 
of this gene resulted increased Mg2+ levels. Increased 
extracellular Mg2+ concentrations caused a decrease in 
the activity of the glutamate Af-methyl-D-aspartate 
(NMDA) receptor [35]. These findings suggest that the 
CNNM2 gene may play an important role in the hypo- 
function of the NMDA receptor, which is implicated in 
the pathophysiology of schizophrenia. 

The inferior frontal gyrus has a multifunctional role in 
human behavior, interpersonal interactions and communi- 
cation [36] . The inferior frontal gyrus consists primarily of 
the heteromodal association neocortex, which is a major 
site of involvement in schizophrenia [37]. Several studies 
have reported that the relative GM in patients with schizo- 
phrenia was significantly reduced in the bilateral inferior 
frontal areas [38-40]. The inferior frontal gyrus can be 
subdivided into three macroanatomical parts: orbital, 
opercular and triangular. The orbital region is one of the 
major regions of the social brain that connects to the orbi- 
tofrontal cortex [41], while the opercular and triangular 
regions form Brocas area, which is an important region 
for speech-language production [42]. We found that the 
CNNM2 genotype affects brain volumes in the orbital re- 
gions of the inferior frontal gyri. Functionally, the orbital 
region is thought to be involved in the processing of em- 
pathy [36] and sentence comprehension [43,44] in the left 
hemisphere and decision-making cognition [41] and fine 
movement control [36] in the right hemisphere. Social 
functions are widely impaired in patients with schizophre- 
nia [45-49]. It is still unclear whether and to what extent 
the effects of CNNM2 polymorphism on GM structure 
observed here might be associated with an increased risk 
for schizophrenia. We suggest that the CNNM2 variant 
may play a role in the social cognition and social function- 
ing impairments noted in patients with schizophrenia 
through GM volumetric vulnerability of the orbital regions 



of the inferior frontal gyri. Further research is needed to in- 
vestigate how a possible relationship between the CNNM2 
gene and hypofunction of the NMDA receptor would result 
in decreased GM volumes of the inferior frontal gyri. 

There was no significant effect of the other four variants 
on any GM volumes. There are several possible reasons 
for the absence of an association. A false negative associ- 
ation cannot be excluded in our study because we applied 
a strict FWE correction for multiple comparisons at the 
whole brain level (p < 0.05). In the Additional file 1: Figure 
S1-S4, the regions shown at the more lenient uncorrected 
threshold of p < 0.001 may be helpful in further studies. 
Interestingly, many of the effects of these genome-wide 
significant variants at the lenient level involved decreased 
GM volumes, including the medial, middle and inferior 
fontal gyri. Reduced GM volumes in these regions have 
been repeatedly demonstrated in imaging studies of 
schizophrenia [14,30]. Another interpretation is that the 
effect of these variants was not sufficiently sensitive to the 
morphological vulnerability of the GM volumes. The effect 
of these variants may be preferable in identifying genotype- 
related vulnerability on other intermediate phenotypes, 
such as cognitive functions and personality traits. There- 
fore, further research is needed to confirm whether the ef- 
fects of these variants could be related to the susceptibility 
of cognitive functioning. 

There were several limitations to this study. Because a 
number of statistical analyses, including the effects of diag- 
nosis, genotypes and their interaction on GM volumes, 
were performed, a correction for multiple testing should 
be considered. However, a consensus on how to correct 
such multiple testing on study combining brain imaging 
and genetics has not been reached in this research field. 
To control type I errors, we applied the strict FWE correc- 
tion for all VBM analyses, while we did not perform any 
correction on the genetic modality. The existence of a false 
positive association cannot be excluded as an explanation 
for our results, although we were quite careful to match 
ethnicity and correct for multiple testing. Further investi- 
gations of other samples with much larger sample sizes 
and/or with different ethnicities and/or in relatives of 
those with schizophrenia are needed to confirm our 
findings. It is unclear whether our results are directly/ 
indirectly linked to the CNNM2 polymorphism rs79 14558, 
to other polymorphisms in high LD with this variant or to 
interactions between this variant and other variants. To 
determine whether rs79 14558 is the most strongly asso- 
ciated variant for schizophrenia and brain structure in 
this gene, an extensive search for other functional vari- 
ants at this locus is needed. Additionally, as with other 
risk variants for schizophrenia, clarifying the biological 
role of this SNP through in vitro and in vivo studies is 
important to improve the understanding of the patho- 
physiology of schizophrenia. 
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Conclusions 

We found that a genome-wide supported variant of 
CNNM2 could be associated with GM morphological vul- 
nerability of the bilateral inferior frontal gyri. These results 
suggest that there may be possible deleterious effects of the 
risk G-allele at CNNM2 in the inferior frontal gyri, which 
may, at least partially, represent the mechanism by which 
CNNM2 increases the risk for schizophrenia. Further re- 
search will be needed to clarify the function of the risk 
CNNM2 variant on the pathophysiology of schizophrenia. 
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